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Abstract 


A method for incorporating the effect of the resonant continuum into 
Hartree-Fock+BCS equations is proposed. The method is applied for the 
case of a neutron-rich nucleus calculated with a Skyrme-type force plus a 
zero-range pairing interaction and the results are compared with Hartree- 
Fock-Bogoliubov calculations. It is shown that the widths of resonant states 
have an important effect on the pairing properties of nuclei close to the drip 
line. 
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In nuclei far from the /3-stability line there are circumstances in which the resonant 
part of the particle continuum plays an important role. One of the first examples was 
given by Migdal |^, who anticipated the existence of halo nuclei like ^^Li by showing that 
two neutrons in a potential well can form a dineutron state in the region of the nuclear 
surface if the potential has a resonant level with energy close to zero. Here, the role of 
the resonance is to enhance the pairing correlations of the two neutrons whose attraction 
is otherwise too weak to bind the system. The interplay between the resonant continuum 
and pairing correlations can also be important for the estimation of unbound processes like 
single-particle excitations lying above the particle emission threshold, which may be found 
especially in nuclei close to the drip line. Thus, if the excitation energy is close to the 
energy of a narrow resonance then one expects an excitation of one-quasiparticle type, in 
which the excited nucleon stays a hnite time before decaying. Besides the unbound processes 
mentioned above, resonant states are important for determining the pairing properties of 
the ground state of bound nuclei far from the /3-stability line. Although in such calculations 
one should consider in principle the complete particle continuum, the largest contribution 
to the pairing correlations is expected to come from the resonant continuum part . 

The aim of this paper is to propose a method for incorporating the effect of the reso¬ 
nant continuum in the Hartree-Fock-|-BCS (HF-I-BCS) approximation. More precisely we 
investigate here the effect of the width of the resonances on the pairing properties of nuclei 
far from stability. As discussed below, this effect is difficult to estimate in self-consistent 
Hartree-Fock-Bogoliubov (HFB) calculations which are presently used for describing pairing 
correlations in nuclei close to the drip line. 

In order to derive the BCS equations in the presence of the continuum we discretize hrst 
the single-particle continuum spectrum by enclosing the nucleus inside a box of a very large 
radius Rf,. This is only a formal step, since will be seen that the parameter Rf, does not 
appear in the hnal results. Thus the genuine continuum is replaced by a set of discrete states 
with the level density given by [^: 
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where gy{e) = {l/TT){2i^ + l){d5u/de) is the so-called continuum level density and 5y is the 
phase shift of angular momentum v = The quantity is the level density in the 

absence of the mean field and is given by = (l/ 7 r)( 2 jy -|- l){dk/de)Rb, where k is the 

momentum corresponding to the energy e. The wave functions corresponding to the positive 
energy states and normalized within the box are dehned by - 01 ,( 6 ,r) = A“^/^(e)</ 3 y(e,r), 
where A0(e) is the norm of the scattering state r) in the box volume. It can be easily 
shown that for the scattering states selected by the box A0(e) = (2jj, -|- 

The gap equations for the states in the box can be written in terms of level density as 
follows: 
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where the indices i,j run over the bound states and is an energy interval associ¬ 
ated with each partial wave The matrix elements of the interaction involv¬ 
ing states in the continuum are given by KeuiyeW = 

{'^p^{e)2plF{e)\V\^jJl{€')'4’-;p■{e')). The rest of the notations are standard p. It can be noticed 
that according to the BCS approach the generalized gap equations above take into account 
only pairing between time-reversed continuum states ue, W. A more general pairing be¬ 
tween continuum states at neighbouring energies is conceivable and this would just be taken 
care of by a continuum HFB approach. At the moment there exist only HFB calculations 
performed with a box boundary condition and we shall numerically compare their results 
with those obtained in the present approach. 

The largest contributions to the integrals above come from the regions where the wave 
functions have a large localization inside the nucleus. This condition is fulhlled in 

energy regions where the F-matrix has poles close to the real energy axis, i.e., near narrow 
single-particle resonances. The integrals can also have large contributions from energy inter¬ 
vals close to zero energy if the S'-matrix has poles corresponding to loosely bound states or 
virtual states near threshold . In the A intervals dehned above the positive energy wave 
functions have the largest localization inside a sphere of radius D, where D is of the order of 
a few times the nuclear radius. Within D the positive energy wave funtions can be related 
to the scattering wave function at resonant energy through simple factorization formulas 
EM- Following Refs. wave function inside D can be approximated by 

r) ^ r) = ry2(e)0^(e^, r) , (4) 


where is the scattering wave function calculated at the resonant energy and nor¬ 

malized within a sphere of radius D. This factorization relation is very useful for evaluating 
matrix elements of hnite range interactions for which it is sufficient to carry space integrals 
over the volume inside the radius D only. For instance, we will use ~ 

and ,1771 1 where on the right hand sides the matrix ele¬ 

ments of the interaction are calculated using the wave functions 4>y{ey,r). For a discussion 
of the accuracy of these approximations see Ref. [§. 

With the help of this factorization the gap equations become: 
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with A,^(e) = ry(e)Aj,. The last expression can be written as g,^{e)A,^{e) = gy{e)Ai, and gives 
the connection between the gaps calculated with the wave functions 'ipp{e) and 4>u{^)- One 
can get the same relation if one writes the gap equation (|^) in terms of a local pairing held 
A(r) of hnite range which cuts oh the contributions of the tail of the wave function 'ipv{e, r) 
beyond the radius D. Thus, 

fRb rD 

A^(e) = / \'ipy{e,r)\^A{r)drK.Ty{e) |0j.(e,r)|^A(r)dr = rj.(e)Aj. . (7) 
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A similar relation can be derived for the positive energy single-particle spectrum. By using 
these relations it can be seen that the gap equations are independent of the box radius. 
As shown above, this is a consequence of the hnite range of the pairing interaction, which 
is sensitive only to the inner part of the resonant continuum wave functions. 

In the BCS approximation the number of particles is hxed consistently with the gap 
equations by counting the particles distributed in the model space in which the pairing 
interaction is effective. Using the same approximations as for deriving the gap equations 
one gets 

N = + Y. [ 9 u{e)vl{e)de . (8) 


This equation together with the gap equations (||,^) are the extended BCS equations for a 
general (hnite range) pairing interaction including the contribution of the resonant contin¬ 
uum. Each resonance is characterized be the quantity A^, which acts as the averaged gap 
of that resonance. 

The above BCS equations are well suited to be specialized to the approximation of 
constant pairing interaction since the resonant states 0i.(ei/,r) and the bound states have 
rather similar localizations inside the nucleus. Then, one just has to replace in Eqs.(§,^ all 
matrix elements by a constant value G. The corresponding BCS equations are the same as 
the ones of Ref. 0. It is worthwhile to point out that in the constant pairing approximation 
as dehned here one preserves the variation of the matrix elements of the pairing interaction 
over the resonance region. Consequently, as seen in Eq.(^), the gap also changes in the 
resonance region and therefore the corresponding pairing held is not constant in the whole 
space. 

We can now extend the above treatment of the resonant continuum to HF-I-BCS. In the 
case of a Skyrme force this is done by including into the nucleon densities the contributions 
of the positive energy states with energies in the selected intervals Iiy and by using the 
factorization relation (^. Thus the resonant continuum contribution to the particle density 
inside the sphere of radius D reads 
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Similar expressions can be derived for the kinetic energy density T{r) and spin density J(r), 
T(r) ^^n2|V^.(r)|2 + ^(n2),|V0.(e.,r)p , (10) 
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where the hrst sum represents the contribution of the bound states. The above densities 
dehne the mean held and the single-particle spectrum. They depend on the occupation 
probabilities and they are calculated iteratively with the BCS equations, as in the usual 
HF-I-BCS calculations |IU| . 


At this point we would like to comment on the relation between the HF-I-BCS equations 
derived above and the HFB approach. The advantage of the HFB approach to treat processes 
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that involve the continuum part of the nuclear spectrum is that the hnite range of the 
pairing field is explicitely taken into account. Therefore, the particle and pairing densities 
acquire automatically a proper asymptotic behaviour ||TT| , p!^ . In order to preserve the same 
behaviour in the HF+BCS limit one should keep the physical condition of a hnite range 
pairing held, as it is also done in HFB calculations in which the pairing held is not calculated 
self-consistently As seen in Eq.(^, a hnite range pairing held implies a cut-oh in the 


tail of the positive energy wave functions. Without this cut-oh the solution of a HF-I-BCS 
calculation with the positive energy states discretized in a box would correspond to a nucleus 
in dynamical equilibrium with a nucleonic gas and not to the nucleus itself [Q. Generally 
the cut-oh radius may be ambiguous but if one restricts oneself to the resonant continuum, 
then there is a rather large region outside the nucleus where the resonant wave functions 
have values close to zero before they start oscillating. In this case the HF-I-BCS results do 
not depend sensitively on the cut-oh radius chosen in such region. 

We now apply the above resonant HF-I-BCS approach to a nucleus far from stability, 
namely ®‘^Ni for which HFB results can be found [0. Here, we wish to compare the results 


of three types of calculations: A) the resonant HF-I-BCS approach where the widths of single¬ 
particle resonances are taken into account; B) a discrete version of it where the widths are 
set to zero; C) the HFB approach of Ref. where the coordinate space equations are 
solved with a box boundary condition and therefore the width ehects would be missing. 
The HF held is calculated with the SHI interaction whereas for the pairing channel 
a zero-range density-dependent interaction is used, H(ri,r 2 ) = Vo(l — p{fi)/pc)5{fi — 7 ^ 2 ), 
where p{f) is the total density, Vq = 1128.75 MeV fm^ and pc = 0.134 fm“^. All calculations 
are carried out up to a distance D= 22.5 fm (see Eq.(4)) but the numerical results discussed 
in the following do not depend sensitively on the precise choice of D in the range (3 — 4)i?, 
where R is the nuclear radius. The resonant states included in the HF-I-BCS calculations, 
with energies smaller than 5 MeV (the energy cut off used in Ref. |T^), together with the 
last bound state are listed in Table 1. The energy (width Fj, ) of a given resonance is 
extracted from the energy where the derivative of the phase shift reaches its maximum (half 
of its maximum) value. The energy intervals R are dehned such that |e — ej,| < 2Ty. 

Let us hrst look at the results of case A. The total averaged gap and the Fermi energy 
are < A >= 0.51 MeV and A=-0.874 MeV, respectively. The total binding energy and 
the pairing energy are £'=-652.7 MeV and £p=3.4 MeV. The corresponding pairing field 
is shown in Fig. 1 while the averaged occupation probabilities and the averaged gaps of 
resonant states and the last bound state 3si/2 are given in Table 1. The change of the 
particle density due to pairing correlations is shown in Fig. 2. It can be seen that in the tail 
region the contribution of the bound states to the total density, given mainly by the loosely 
bound state 3si/2, is dominant. In order to see how the neutrons are distributed at large 
distances, we have calculated the number of neutrons outside a sphere of radius 12 fm. We 
hnd that the total numbers of neutrons distributed in bound and resonant states between 
12 and 22 fm are 0.069 and 0.037, respectively. 

A proper estimation of the particle distribution at large distances is difficult in HFB 
calculations based on a box boundary condition. Due to the box the wave functions which 
are spread far from the nucleus are generally pushed towards smaller distances. Thus, in the 
HFB calculations of Ref. [Q the occupancy of the loosely bound state 3si/2 which gives the 
dominant contribution in the density tail region, depends strongly on the box radius and is 
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generally underestimated. Consequentely the HFB density is smaller in the tail region than 
the HF+BCS density. As seen in Fig. 2 the tail of the HFB density is actually smaller even 
than the HF density. 

For case B the total averaged gap becomes < A >= 0.72 MeV and the Fermi energy is 
A= -0.948 MeV. The binding energy increases to the value E=-653.1 MeV while the pairing 
energy is equal to Ep=Q.2 MeV. From Table 1 one can see that the occupancy of the resonant 
states is almost doubled compared to the case when the effect of the widths is taken into 
account. The corresponding changes in the pairing held and particle density are shown in 
Figs. (1,2). As seen in Fig. 2 the pairing held given by the HFB calculations based on a box 
boundary condition are quite similar to the HF-I-BCS pairing held of case B, i.e., calculated 
without taking into account the width ehect. The same similarities are seen in the Fermi, 
binding and pairing energies, which in the case of HFB are: A=-0.956 MeV, i7=-653.7 MeV 
and Ep=6.9 MeV. These similarities show that in the HFB calculations with a box boundary 
condition the resonant continuum is actually described by quasibound states and therefore 
in such calculations the ehect of the widths of resonant states upon pairing properties is not 
taken into account properly. 

In conclusion one sees that neglecting the contribution of the widths in HF-I-BCS cal¬ 
culations one enhances artihcially the amount of pairing correlations. This enhancement is 
due to the fact that if the width is neglected then all the pairing strength is collected from 
the scattering state at the resonance energy. At this energy the scattering wave function has 
the highest spatial concentration (as compared with the nearby scattering functions) within 
the nuclear region. The effect of the width is to diminish the pairing strength because the 
pairs can now scatter also in the nearby states around the resonance energy which are less 
conhned inside the nucleus. In a time dependent picture the dependence on the width seen 
above would correspond to the fact that a pair scattered to a resonance state has a hnite 
lifetime, thus contributing less to the pairing correlations as compared to the case of a very 
narrow (quasibound) state. In the HF-I-BCS approach presented here this effect is taken 
into account automatically through the continuum level density. 

Formally, in a coordinate space HFB approach the contribution of the whole continuum 
is taken into account and therefore the effect of the resonant continuum discussed above 
should be also present. Furthermore, a HFB calculation contains also correlations from 
pairs in states which are not time-reversal partners. These correlations which are absent 
in a HF-I-BCS approach could be important for the scattering states with energies close to 
a resonance because their wave functions have similar localization properties [^. However, 
the estimation of the effect of the resonance widths upon pairing correlations is still an open 
problem in the existing self-consistent HFB calculations. This is because the numerical 
methods used for solving the coordinate space HFB equations are based on discretizations 
of the continuous spectrum. With the currently used values of box radius (R ~ 15-30 
fm) each resonance is represented by a single discrete state in the spherically symmetric 
case. To have a discrete level density high enough to describe properly the shape of the 
resonance would require an extremely large box radius, a condition which makes practical 
calculations untractable. In addition, the use of a box of hnite radius could alter signihcantly 
the description of the pairing correlations induced by the loosely bound states whose tail 
might extend beyond the box radius. Although the coordinate space HFB approach is in 
principle the appropriate tool for treating continuum effects, solving the HFB equations by 
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imposing a box boundary condition to the solutions does not guarantee that all effects of 
the continuum, particulary the effect of the widths of resonant states, are properly taken 
into account. 

In summary, a method to include the resonant continuum in the HF+BCS approximation 
is presented. We have here concentrated on the regions of the continuous spectra which are 
close to single-particle resonances because they bring the most important contributions to 
pairing correlations. The method can be used to take also into account the effect of non¬ 
resonant continuum states close to the continuum threshold, which can be important in the 
presence of loosely bound states or virtual states. In the numerical example it has been 
shown that the widths of resonant states have an important effect on the pairing properties 
of nuclei close to the drip line. In order to describe this effect in a self-consistent HFB 
approach one needs to solve the HFB equations with proper boundary conditions. This 
work is in progress. 

We thank P.H.Heenen and J. Terasaki for useful discussions and for the detailed results 
of their HFB calculations. Two of us (NVG and NS) would like to thank the IN2P3-INPE 
collaboration in the frame of which part of this work was performed. NS acknowledges the 
hnancial support of the Wenner-Gren Foundation. 
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TABLES 


TABLE I. Results of HF+BCS calculations for the nucleus ®^Ni. and are the averaged 
gap and averaged occupation probability of the single-particle state n of energy e„, and width 
En- The notations A^, v^, Sn and stand for the corresponding quantities calculated without 
including the effect of the widths of resonant states in HF-I-BCS equations. The single-particle 
energies, their widths and the pairing gaps are expressed in MeV. 


n 







A, 

An 

Sl/2 

-0.647 

-0.644 

— 

— 

0.295 

0.294 

0.505 

0.674 

ds/2 

0.441 

0.417 

0.077 

0.068 

0.041 

0.075 

0.630 

0.847 

9712 

1.604 

1.587 

0.009 

0.008 

0.029 

0.055 

0.966 

1.306 

hii/2 

3.309 

3.302 

0.017 

0.016 

0.017 

0.034 

1.227 

1.658 


FIGURE CAPTIONS 


Figure 1: Neutron pairing field as a function of the radius. The full (dashed) line shows 
the results of the HF-I-BCS calculations with (without) the effect of the width included. The 
line marked by crosses shows the HFB results of Ref. 0. 

Figure 2: Neutron density in ®^Ni, calculated in HF (long dashed line) and HF-I-BCS. 
In the order of decreasing tail the results of HF-I-BCS correspond to the following approx¬ 
imations: effect of the widths neglected; effect of the widths included; contribution of the 
bound states to the density. All densities are in fm“^. The inset shows the corresponding 
densities in linear scale. 
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